Optimally Controlled Field-Free Orientation of the Kicked Molecule 
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Efficient and long-lived field-free molecular orientation is achieved using only two kicks appro- 
priately delayed in time. The understanding of the mechanism rests upon a molecular target state 
providing the best efficiency versus persistence compromise. An optimal control scheme is referred 
to for fixing the free parameters (amplitudes and the time delay between them). The limited number 
of kicks, the robustness and the transposability to different molecular systems advocate in favor of 
the process, when considering its experimental feasibility. 
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Molecular orientation, in particular during the free 
evolution of the system, has recently been abundantly 
discussed in the literature as a process playing an im- 
portant part in a variety of laser-molecule control issues, 
amoiig which are chemical reactivity 0], nanoscale de- 
sign 0, 13 , surface processing 0, |3 , and attosecond time 
scale pulse production 0, Q. A basic mechanism by 
which orientation is achieved involves sudden optical ex- 
citation, such as half-cycle pulses (HCP), that impart a 
kick to the molecule which orients itself along the polar- 
ization vector of the linearly-polarized electromagnetic 
field (the kick mechanism) _10] . 

Application of a series of kicks is in fact a gen- 
eral control strategy that can enhance molecular align- 
ment 01 or even lead to a squeezing of atoms in an 
optical lattice IT^. Moreover, alignment by a pair 
of pulses has been experimentally achieved ^a, [ig • The 
generic strategy, as first suggested in Ref. rests on 
the application of a sudden impulse when a certain ob- 
servable reaches a maximum (or minimum), such as the 
expectation value of cos 9 (which measures orientation, 9 
being the angle positioning the molecular axis with re- 
spect to the laser polarization) for the control of orienta- 
tion. 

In this paper we present a different approach based on 
a target state 18j, instead of an observable, which al- 
lows to consider not only the efficiency of the orientation 
((cos 6*)), but also its persistence. This stems from the 
fact that, when dealing with field- free orientation of a 
molecule, a compromise between efficiency and duration 
has to be looked for in the optimization criterion 1^. 
Although the generic strategy (i.e., kicking when an ob- 
servable reaches its maximum value) can be applied also 
to reach the target state 20j, we show here, that this 
strategy is far from optimal, and that two kicks can be 
sufficient to come within 1% of the target state. 



A. Model 

A diatomic molecule illuminated by a moderate ampli- 
tude HCP is described within a model of a rigid rotor in- 
teracting with the field through its permanent dipole (po- 
larizability interaction neglected). The time-dependent 
Schrodinger equation (TDSE) accounting for the time 
evolution of the system is (in atomic units) 



BJ^-flQ-e£{t) i^{9,ip;t), (1) 



where the molecule is described by its rotational constant 
B and permanent dipole moment /Iq, and the HCP is 
characterized by its amplitude £{t) and unit polarization 
vector e*. J is the angular momentum operator. The wave 
function involves two angular variables: the polar 9 and 
azimuthal ip angles. Due to cylindrical symmetry around 
e, the motion associated with the azimuthal angle ip is 
completely separated and is not considered hereafter. As 
a consequence, the projection quantum number M of J 
on e'is fixed. The polar angle 9 takes part (apart from the 
analytical characterization of J^) in the dot product fig -e 
in Eq. through its cosine. The short duration r of 
the HCP, as compared to the molecular rotational period 
Trot — T^/B, is accounted for through a dimensionless 
small parameter e = T/T^ot — tB/tt less than one. A 
time scaling is achieved ^21, 22] by introducing a variable 
s = t/Tiot (s G [0,e] during the HCP pulse), leading to 
the following form of the TDSE iPl 



i—'4){9,if\s) 
as 



ttP _ E{s)cos9 ilj{9,ip;s) 



(2) 



where E{s) — 7r^£(Ti-ots)- A sudden approximation, 
where the small perturbation parameter is the duration 
of the pulse e, provides the wave packet at time s, as a 
result of successive applications of two unitary evolution 
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operators on the initial state tp{9, (p;s = Q) (taken as a 
pure quantum state), 



ij{0,ip;s)^Uois)U^iJj{0,ip;s = O), 



where 



and 



Uq{s) = exp 



Ua = exp [iA cos 9] 



(3) 
(4) 
(5) 



with A = E{s)ds a dimensionless parameter combin- 
ing molecule and field characteristics, integrated over the 
whole pulse duration. The initial molecular pure quan- 
tum state is taken as the isotropic spherical harmonic 
Yj^M=oi0,v)- The dynamical picture, which emerges 
from Eqs. for a single HCP, is that a unitary oper- 

ator Ub. imparts a kick, measured in terms of a strength A 
times cos 9, to the molecule in its initial state. This pro- 
duces a rotational excitation bringing the molecule from a 
completely isotropic angular distribution to a more ori- 
ented configuration. The further field-free evolution is 
monitored by Uo{s). In the case where a series of HCPs 
is considered, the evolution operators UqIs)!/^ have to be 
applied for each individual pulse, taking into account the 
time delays between each. It is also worthwhile noting 
that introducing the dimensionless parameters e, A and 
variable s helps to free oneself from a specific molecule- 
plus-field system. Actually, molecular and field charac- 
teristics are combined in such a way that a large rota- 
tional constant B can be overcome by a shorter pulse 
duration r, the relevant condition being e <C 1. Simi- 
larly, a small permanent dipole /io can be overcome by 
a stronger field amplitude £, the relevant parameter A 
being proportional to the product fio£/B. As for the 
rescaling of time, it gives access to results in terms of a 
dimensionless time s that is taken as a fraction of the 
rotational period, which again is molecule independent. 



B. Target state 

A measure of orientation is given by the dynamical 
expectation value of cos0, i.e., 

(cos 6*) (s) = {'iP{9, ip; s)\ cos9\^j{9, ip; s)) , (6) 

the optimization of the orientation being in relation with 
the maximization (or minimization) of {cos9) (s). We 
have recently developed a generic strategy that, when 
applied to cos9 taken as an operator, can be summarized 
as following [l8l |: 

(i) The physical Hilbert space in which cos 9 is acting 
is reduced to a finite subspace Ti^^'^ of dimension N. The 
expectation value of the projection of cos^ on this sub- 
space is time periodic and can be represented, in the basis 
of spherical harmonics, as a finite N x TV-dimensional ma- 
trix with discretized, bounded eigenvectors (as opposed 
to cos 9 which has a continuous spectrum) . 



(ii) A target state can be defined as the eigenstate of 
the projection of cos 9 on Ti^^) corresponding to its high- 
est (or lowest) eigenvalue. Such a state can be explicitly 
calculated, in the basis of spherical harmonics, by diag- 
onalizing the corresponding matrix. The full advantage 
of the dimensionality reduction remains in the fact that 
a smaller N involves a lower rotational excitation that 
allows for a longer duration of the orientation after the 
pulse is over (see Fig. 1 of Ref. 0|). The most excit- 
ing observation is that a target state calculated within 
a subspace as small (with respect to dimensionality) as 
N = 5 already leads to an excellent orientation efficiency 
of {cos 9) ~ 0.9. In other words, the target state fulfills 
the two requirements of the orientation control problem, 
and in that respect is far superior to all other intuitive 
criteria that have been previously used |l9l | . 

(iii) A strategy is proposed to reach this target state 
by applying a series of identical short pulses at times 
when (cos 9) reaches its maximum (or minimum) follow- 
ing field-free evolution. Consequently, the corresponding 
wave function converges to the target state. The robust- 
ness of the strategy has been checked against the pulse 
strength and the time delays. For completeness, we have 
also to mention a similar strategy (leading to similar re- 
sults) that consists in applying the pulses every time the 
projection of the time-evolved wave packet reaches its 
maximal projection on the target state. 

From this theory, that can actually be generically 
transposed to other control issues, the recipes that 
emerge for a possible control of molecular orientation is 
the application of a train of short and identical pulses 
with a given total area (within 10 to 15% of accuracy) and 
respecting predetermined time delays between successive 
pulses (within 10% of accuracy). For typical cases, 10 to 
20 pulses with A = 1 are necessary to reach the target 
state. But obviously this strategy taken as a whole (with 
the values of the time delays, the number of pulses and 
their integrated amplitude A) is not unique. However, 
it indicates that a train of short, time-delayed pulses 
induces a repeated kick mechanism at specific molecu- 
lar response times that improves the efficiency and the 
duration of orientation. It is precisely this information 
that serves here as a basis for a numerical optimal con- 
trol scheme conducted using an evolutionary strategy 
(ES) H. 

The target state being clearly identified, the optimiza- 
tion aims to maximize the projection of the instanta- 
neous wave packet on the target. More precisely, the 
wave packet at time s being expanded on the basis of 
spherical harmonics yj,M=o(^, </') as 



V^(0,^;s) = ^cj(s)y7,M=o(^,^), (7) 

J=0 

the probability to be maximized is 

P(,s) = |c*(s).Ctargct|', (8) 

where c*(s) = (cq(s), c|(s), . . . , ) and Ctargot is the cor- 
responding column vector of the weighting coefficients of 
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the target state on the same basis of spherical harmonics. 
The ES deal as usual with the minimization of a criterion 
defined here as 

J = l-P(s,(P^ax)), (9) 

where Si(Pniax) is the time for which P reaches a max- 
imum during the free evolution over a rotational period 
following the radiative interaction. The parameters of 
the optimization are the amplitudes {A) and time inter- 
vals characterizing a train of kicks that can be produced 
by HCPs. The ES is implemented using the Evolving 
Objects library (EOlib) ,ii,,2^. 



III. RESULTS 

All the calculations that are presented here deal with 
a target state in a A'^ = 5-dimensional finite subspace 
H^^^{the maximum allowed rotational excitation J,nax = 
4 being the one previously retained as satisfying the 
best post-pulse orientation efhciency/duration compro- 
mise ilSi]). The wave packet by itself is propagated in 
a larger, although finite, subspace and thus may reach 
higher rotational states. We note that the use of a tar- 
get state defined in a reduced Hilbert subspace dispenses 
from having to penalize with respect to the total pulse 
intensity. Too strong pulses will necessarily move the 
system outside the subspace and thus reduce the value 
of the projection on the target state. Nevertheless, us- 
ing a target state is not equivalent to penalizing since, 
for the same pulse energy, there might exist states that 
show better orientation but with less persistence. Only 
the target state ensures that the efficiency/persistence 
compromise is achieved. The optimal control strategy is 
guided by two different approaches, depending upon the 
parameter space chosen. 



A. Time delays as the only parameters 

For a given number of pulses, the only task conferred 
to the ES algorithm is the optimal determination of the 
time delays between the pulses for a minimization of j 
[Eq. 0]. The pulses are considered identical, with all 
other parameters taken such that A = 1. The results 
for 3 and 4 pulses are displayed in tables U and m col- 
lecting the values of the time delays Si between pulses i 
and i + 1 and their times of application Si , together with 
the values reached for P(si(Pinax)) after each successive 
kick. Are also given the time intervals (5i(coSmax) and 
Si{Pmax) between the ith kick and the next maxima of 
{cos 9) (s) and of P(s). respectively, for comparison with 
the strategies of Ref. [T^. Concerning the three-pulse 
model of Tab. two observations can be made, show- 
ing that the general theory of Ref. summarized 
in Sec. IIIBI is neither unique nor optimal. A value of 
P = 0.8905 is reached using the optimal time delays Si, 



i 




Si 


(5i(cOSmax) 


(-fmax ) 


p 

max 


1 


0.1065 


0. 


0.2070 


0.1607 


0.3913 


2 


0.0029 


0.1065 


0.0991 


0.0855 


0.6685 


3 




0.1094 






0.8905 



TABLE I: Optimal delays 5i (and corresponding times be- 
tween Si) between 3 kicks of constant amplitude A = 1, along 
with the projection Pmax on the target state. (5i(coSmax) and 
^^(fmax) are the time intervals between the ith kick and the 
next maxima of (cosS) (s) and of P(s), respectively. 



i 


5^ 


Si 


5i(cOSmax) 


(-frnax) 


p 

max 


1 


0.2388 


0. 


0.2070 


0.1607 


0.3913 


2 


0.9901 


0.2388 


0.0907 


0.0962 


0.4848 


3 


1. X 10"^^ 


1.2289 


0.0884 


0.0853 


0.8018 


4 




1.2289 






0.9787 



TABLE II: Same as for Tab.Q but with four identical pulses. 



far better than the one that can be reached by applying 
the kicks precisely at the maxima of {cos 6) (s), yielding 
P = 0.5277. The comparison between Si's and (5i(coSmax) 
and (Si (Pmax) 's shows that the optimal strategy is to ap- 
ply the pulses before the maxima of (cos 9) (s) or of P(s). 
The situation is different for the 4-pulse model, Tab. ^1 
advocating again for the non-uniqueness of the solution. 
The second pulse in particular is applied after a time 
delay Si larger than (5i(coSniax) or (5i(Pmax)- An excel- 
lent value for P is obtained, showing that after 4 pulses 
the molecular state that is reached is close to the tar- 
get within 2%. But, even more interestingly, the third 
and fourth pulses are applied at times very close to one 
rotational period (i.e., S2 = 0.99, 5^ = 0.00) after the 
second kick. Due to the periodicity, this amounts to ap- 
plying simultaneously three identical pulses after a time 
delay Si corresponding to the second pulse. Still another 
way of analyzing the situation consists in applying a first 
pulse with an amplitude Ai = 1.0 and after a time de- 
lay of Si = 0.2388 applying a second pulse of amplitude 
A2 = 3.0. Such a strategy has actually been checked and 
leads to a final projection P = 0.9776, very close to the 
one displayed in Tab.|lTl i.e., P = 0.9787. 

We note that a similar result is also obtained with five 
kicks, allowing then to reach P = 0.9928. Better results 
are not reached with more kicks as the subsequent ones 
have the counter effect of increasing the rotational excita- 
tion, and thus of pushing the system outside the Hilbert 
space 7i'^' where the target resides. 
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order to reach a target state corresponding to an oriented 



I I I I I I I I I I I I I I I I I I I I I I I 
0.0 0.2 0.4 0.6 0.8 1.0 1.2 

FIG. 1: Resulting orientation, as expressed by {cos 9), for 
the optimal 2-kick solution, along with the amplitudes of the 
kicks. 




B. Time delays and amplitudes as parameters 

Referring again to the strategy of implementing in the 
ES what is learned from previous attempts, we extend 
the parameter space such as to account for the variation 
of both the time delays Si and the amplitudes Ai. In 
addition, we restrict the optimization scheme to a 2-pulse 
model, which results into the rather simplified task for 
the ES of providing merely with 3 parameters: Ai, A2, 
and ^1. 

The optimization then yields the strategy of giving a 
first kick of amplitude Ai = 0.9741, followed by the sec- 
ond of amplitude A2 = 3.2930 after a delay Si = 0.2419, 
allowing to reach P = 0.9886. This confirms the previ- 
ous observations: the target state is easily reached us- 
ing only two pulses and the intensity of the kicks must 
be restrained so that highly excited rotational states are 
not populated. The resulting time evolution of (cos 9) is 
given in Fig. ^ A value of (cos 6') — 0.9078 is reached, 
slightly greater than that of the target (0.9062). The cor- 
responding angular distribution at maximum orientation 
is shown in Fig. 13 along with the target state: the two 
are virtually indistinguishable at this scale. 

We have checked the robustness of this strategy by 
varying the parameters by 10% of their optimal value. 
The smallest value of P obtained is then 0.9689, so the 
results remains quite close to the target. Orientation 
efficiency remains within 2.2%, whereas its duration is 
shorten by at most 9%. 



FIG. 2: Polar plot of the angular distribution reached at the 
maximum of orientation, t — 0.3042Trot (solid line) compared 
with that of the target state (dotted line). 



molecule. The solution turns out to be very efficient, al- 
lowing to reach the target state (within 1%) with only 
two pulses, instead of the approximately 15 pulses of the 
previous mathematically built (but not unique) strategy 
ofRef. El for the same system. This advocates for great 
experimental feasibility and, even more importantly, it 
points out the broad interest of the overall methodology. 
As has been already shown, the mathematically clear de- 
piction of a quantum target state in a finite dimensional 
subspace can be conducted for a large class of observ- 
ables (some examples in comparing systems interacting 
with a thermal baths or other dissipative environments 
are provided in Ref. jl^). Once the target state is de- 
fined, ES can be successfully run with a simple criterion 
of maximum projection on the target state referring as a 
basic mechanism to a train of kicks. The only parameters 
to be optimized are the time delays and amplitudes. A 
small number of such kicks, reachable within modest ex- 
perimental environment, allow for a remarkably efficient 
and persistent control. 



IV. CONCLUSION 
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